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We investigate in detail the low-energy spectrum of the R-parity violating minimal supergravity 
model using the computer program SDFTSUSY. We impose the experimental constraints from the 
measurement of the anomalous magnetic moment of the muon, {g — 2)^, the decay b ^ sy as 
well as the mass bounds from direct searches at colliders, in particular on the Higgs boson and 
the lightest chargino. We also include a new calculation for the R parity violating contribution 
to Br(_Bs — > ij,'^fi~). We then focus on cases where the lightest neutralino is not the lightest 
supersymmetric particle (LSP). In this region of parameter space either the lightest scalar tau 
(stau) or the scalar tau neutrino (tau sneutrino) is the LSP. We suggest four benchmark points with 
typical spectra and novel collider signatures for detailed phenomenological analysis and simulation 
by the LHC collaborations. 



I. INTRODUCTION 

The most widely discussed solution to the gauge hi- 
erarchy problem ^ of the Standard Model of parti- 
cle physics (SM) [2[is the supersymmetric extension 
of the SM (SSM) Q. However, no supersymmetric 
particle has been observed to-date. Thus if it exists 
supersymmetry (SUSY) must be broken and the mass 
scale should be of order the TeV energy scale, in or- 
der to maintain the solution to the hierarchy prob- 
lem Q- The TeV energy realm will be probed at the 
LHC starting in 2007 with the search for SUSY of 
paramount interest Q. 

The general SSM renormalisable superpotential with 
minimal particle content is given by 
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Here we have employed the notation of Ref. We 
have split the superpotential into two sets of interac- 
tions, for reasons which we shall explain shortly. 

If simultaneously present, the baryon and lepton num- 
ber violating interactions in Eq. ([3]) lead to rapid pro- 
ton decay Therefore the SSM must be augmented 
by a symmetry. The most widely studied scenario is 
R-parity 0, for which — 0. However, it suffers 
from dangerous dimension-five proton-decay operators 
0. This is solved by proton hexality, Pg, an anomaly- 
free Zg discrete gauge symmetry [10[ . The renormalis- 
able superpotential terms are equivalent to conserved 
R-parity, namely Wpg , but the dimension-five proton- 
decay operators are forbidden. So we consider this the 
more appropriate choice and have correspondingly la- 
belled the superpotential in Eqs. pl^ . The Pg-SSM 
is conventionally denoted the minimal supersymmetric 
Standard Model (MSSM). 

Baryon triality, is an anomaly-free Z3 discrete 
gauge symmetry which prohibits the UDD in- 
teractions, thereby stabilising the proton. It is thus 
theoretically equally well motivated to Pg. A third 
possibility to stabilise the proton is lepton parity , 
which prohibits the lepton-number violating interac- 
tions in Eq. ([3]). However, lepton parity is anomalous 
[ut - As the alternative to the Pg-MSSM we shall focus 
on the phenomenological analysis of baryon triality at 
colliders. However, we shall also consider some impli- 
cations of lepton parity. 

The Pg-MSSM mass spectrum and couphngs have been 
studied in great detail 0, However, there 

are two reasons why it is highly desirable to perform 
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analyses of the general SSM at colliders. First, to 
make sure that no possibility to discover SUSY is left 
untried. If we focus our endeavours on one specific 
form of SUSY, due to theoretical prejudice for exam- 
ple, then we might initially entirely miss the discov- 
ery. Second, SUSY is not a simple extension of the 
SM, like for example right-handed neutrinos, which 
does not affect most parts of the model. In the SSM, 
the particle spectrum is more than doubled. If the 
SSM is to provide a solution to the technical hierar- 
chy problem, we expect all of the new particles to have 
masses ~ 0{2 TeV). Thus if low-energy SUSY exists 
we should expect many new processes to be kinemati- 
cally accessible at the LHC. When considering a spe- 
cific SUSY process as a signature, many other SUSY 
processes will act as background and must be taken 
into account 

Thus in order to perform simulations of signatures 
involving the complete supersymmetric spectrum in 
preparation for the LHC, we restrict ourselves to 
specific points of the Pg-minimal supergravity (Pg- 
mSUGRA) parameter space. Over the years several 
such example points have been proposed [l^ [H, . 
Some of these "benchmark points" have been excluded 
by searches at LEP, for example the points 3 and 
5 in Ref. [Tsj . Some points are strongly motivated 
by the resulting cosmological relic density [ll]. A 
set of benchmark points and parameter lines in the 
mSUGRA parameter space have been agreed upon, 
the Snowmass Points and Slopes (SPS) The re- 

sulting spectra have been compared for various nu- 
merical RGE integration codes jl^, US] and have been 
investigated both in the context of the LHC and the 
ILC. It is our purpose to facilitate the extension of 
this work to the Pe-MSSM. As a first step, we shall 
investigate the mass spectrum and couplings in some 
detail. For further more detailed experimental studies 
we also propose a set of benchmark points. 

In Sect. [TT] we briefly review the most widely studied 
model, namely the Pg-MSSM embedded in supergrav- 
ity [U, [13 • In order to throw a wider net in our 
search for SUSY, we extend this discussion to the Jg- 
MSSM in Sect. IIII|, emphasising the main differences 
to the Pg conserving case and examining bounds on 
the couplings. In this first study we sometimes 
will focus on the simplest model of no-scale super- 
gravity [2^. We examine constraints on ^g mSUGRA 
parameter space arising from precision measurements 
in section HVl and the sparticle spectrum in section fVl 
We present the definition of the proposed mSUGRA 
benchmark points in section IVI) where we also de- 
tail the spectrum and sparticle decays. We summarise 
and conclude in section IVIIl following with expres- 
sions for the ^g contribution to Bi-{Bq^ e+Cj") de- 
cays, a generalisation of results in the literature in Ap- 
pendix]^ In this paper we do not consider alternative 
mechanisms to communicate to the visible sector such 
as gauge or gaugino mediated SUSY breaking [24] or 



anomaly mediated SUSY breaking [2£ 



II. THE Pb-MSSM 

The Pg-MSSM has 124 free parameters [11], and the 
|*g-MSSM more than 200. Such an extensive param- 
eter space is intractable for a systematic phenomeno- 
logical analysis at colliders. It is thus mandatory to 
consider simpler models, which represent the variety 
in the SSM. Hereby both observational hints [l^] and 
theoretical considerations '28] have been employed as 
guiding tools. 

Within the SSM, the unification of the three SM gauge 
couphngs at a scale Mx = 0(10^'' GeV) [27] indicates 
the embedding in a unified model. The most widely 
studied and the sirnplest of these models is mSUGRA 
with conserved Pg 281 ] and radiative electroweak sym- 
metry breaking [29). In this case, we are left with five 
free parameters 



, Ml/2 , ^0 , tan/3 , sgn(^) 



(4) 



where the first three impose boundary conditions on 
the soft-SUSY breaking parameters at Mx, while the 
other two are boundary conditions at Mz- Here Mq is 
the universal, i.e. flavour independent, soft breaking 
scalar mass. We thus have for example 



Mo, @Mx , (5) 



where m,~,, denotes the sneutrino mass, denote 
the 5[/(2)-doublet (L) and SU{2) singlet (R) charged 
slepton masses, and m^j^ ^ denote the squark masses. 
™-ffi.2 are the soft breaking Higgs boson masses. M1/2 
denotes the universal gaugino soft breaking mass and 
we have for the bino (Mi), the SU (2)-wino {M2) and 
the gluino mass (M3), respectively [30| 



Ml = M2 = M3 = Ml/2, 



@M 



X 



(6) 



Aq is the soft breaking universal trilinear scalar in- 
teraction [3l| . tan f3 — ^ is the ratio of the vacuum 
expectation values of the two Higgs doublets Hi 2 and 
sgn(/i) is the sign of the Higgs mixing parameter of 
Eq. ([2]). The absolute value of /x as well as a possible 
bilinear scalar interaction Bq are fixed by the elec- 
troweak symmetry breaking conditions at the SUSY 
breaking scale [29| . 

Given the five parameters in Eq. fl]), we can com- 
pute the full supersymmetric spectrum and couplings 
at the electroweak scale through the renormalisa- 
tion group equations (RGEs). This has been stud- 
ied in great detail [H, [H, [l3]. An example spec- 
trum is shown in Fig. [1] for the SPSla parame- 
ter set [Mo = 100 GeV, M1/2 = 250 Ge)^ = 
-100 GeV, tan/3 = 10,sgn(/x) = +1] In 
the left column (red), we have the neutral scalars 
{h^ , , as well as the charged Higgs bosons 
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FIG. 1: The SPSla supersymmetric spectrum. The Pe- 
mSUGRA parameter values are Mo = 100 GeV , M1/2 = 
250 GeV,ylo = -100 GeV , tan /3 = 10,sgn(^) = +1 . 
The up, down, charm and strange squarks are quasi- 
degenerate in mass and are labelled by qL,R- 



(H^). In the second column (black), we have 
the charged slepton masses. In the third column 
(light blue), we show the neutrahnos (x?=i,2,3,4) ^^'^ 
charginos (xf2) ^^'^ the fourth column (marine 
blue) the squarks (w, d, c, s, t, b) and gluinos g. 

We next discuss some typical features of this spectrum. 
First, the lightest supersymmetric particle (LSP) is 
the lightest neutralino. The SPSla parameter values 
have been specifically chosen to ensure this. Since the 
LSP is stable in the Pg-MSSM, it must be electrically 
and colour neutral for cosmological reasons [12] , and a 
stable sneutrino LSP is experimentally excluded [ssj . 
In Pe-mSUGRA, the regions of parameter space where 
the Xi is not the LSP are thus excluded from further 
consideration. This is important for what follows since 
in J*g-mSUGRA this region is re-opened. In Fig. ^ we 
have plotted the nature of the LSP in the Mq-Mi/2 
plane around the SPSla point. Most of the xO— LSP 
parameter space leads to dark matter relic densities 
that are too high compared with the WMAP3 con- 
straints if Pg is conserved. Since we will consider the 
^5 case where the LSP decays and does not constitute 
dark matter, we do not plot the constraint on Fig. [21 
The blackened out region at smaller values of Mi/ 2 is 
excluded due to the LEP2 Higgs exclusion bound of 
m/j > 114.4 GeV at 95%C.L. (3^], and/or the require- 
ment of the absence of tachyons in order to obtain a 
valid electroweak vacuum. Note that, anticipating a 
3 GeV error in SOFTSUSY's prediction of mh, we have 
imposed a lower bound of 111.4 GeV in the figure. In 
large regions of parameter space the Xi is indeed the 
LSP. Separated from this by the black contour, for low 
values of Mq and large values of M1/2 the lightest f is 
the LSP. For later use, we also display the stau mass 
in the figure. 



Second, in Fig. [T] it can be seen that the squarks and 
gluinos have masses much larger than Mq and are 
much heavier than the sneutrinos and charged slep- 
tons. This is because the strong interaction dominates 
the RGE running of the former. For Mq :s> M1/2 the 
slepton masses can be of order the squark and gluino 
masses, cf. SPS2 [17]. Third, one stop, the one which 
is predominantly in, is significantly lighter than the 
other squarks since the stop mixing sin 6t is of order 
TUt/MsuSY, whereas Msusy is the SUSY breaking 
scale. Similarly the stau, which is dominantly f^, is 
the lightest charged slepton, although the relative ef- 
fect is much smaller, since rrir ^ Mq. Since the slep- 
tons themselves are relatively light the ffl is the next- 
to-lightest supersymmetric particle (NLSP). Fourth, 
the lightest Higgs boson is significantly lighter than 
the other Higgs bosons. 

For SPSla and other typical mSUGRA points the 
largest supersymmetric production cross sections at 
the LHC are for gluinos and squarks ^1]. Usually, 
squarks and gluinos cascade decay via intermediate 
SUSY states to the LSP, leading eventually to a large 
missing pT signature of the Pg-SSM. 

In order to prepare for the LHC, it is important, given 
a supersymmetric spectrum, to analyse the resulting 
production processes and decays in great detail. In 
particular, this involves the full detector simulation for 
each signature. Even with the highly reduced parame- 
ter space of Pg-mSUGRA, this is beyond present capa- 
bilities. There is a special mSUGRA model with even 
fewer parameters: the so-called no-scale supergravity 
model 1371 , where we have as additional conditions 



Aq = Mq = , 



X ■ 



(7) 



This is however experimentally excluded through a 
combination of the LEP Higgs and chargino search [s^] 
and the requirement of a neutral LSP. On the other 
hand, in the ^g-SSM the LSP need not be neutral 
and regions of no-scale mSUGRA parameter space are 
allowed [33. 



III. THE Pe-VIOLATING MSSM 

We next investigate the mass spectrum of the f^- 
MSSM with particular focus on phenomcnological 
analyses in preparation for the LHC. The superpo- 
tential is given in Eqs. (H])-® and has 48 additional 
parameters, Xijk, ^'ijk-i ^'Ijk-i '^i' with further param- 
eters when we include SUSY breaking. Again, this 
parameter space is too large for a detailed phenomcno- 
logical analysis and we restrict ourselves to the sim- 
plification of minimal supergravity with universal soft 
SUSY breaking parameters at Mx- The ^-mSUGRA 
model was developed in Refs. [H, [H, [s^, E^. The 
RGEs connecting Mx with the weak scale are given at 
one-loop in Refs. [H, IH, [13, El . The two-loop RGEs 
for the SM gauge couplings and the superpotential pa- 
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FIG. 2: Parameter space around SPSla. The bar on the 
right displays the lightest stau mass. Aq — —100 GeV, 
tan/3 = 10, sgn{fi) = +1. The blackened out region to the 
left is excluded due to the LEP bounds and tachyons. The 
black contour distinguishes between areas with x^-LSP and 
f-LSP. 

rameters [i^l can have a qualitative impact on gauge 
coupling and b-r unification and are included in our 
analysis. The impact of three-loop RGE terms is usu- 
ally 1-2% [i^l for the squarks, and negligible for weakly 
interacting sparticles. The two- loop RGEs for the soft- 
breaking terms have recently been calculated [1^ . We 
restrict ourselves to one-loop RGEs for the soft break- 
ing parameters in order to perform scans in parame- 
ter space in a reasonably small amount of CPU-time. 
This is likely to be an excellent approximation for fg 

couplings -^0.1. 

In order to keep the parameter space manageable, we 
shall also only consider one non-zero parameter of Wy^ 
at a time at Mx- However, we include the possible 
dynamical generation of the other parameters through 
the coupled RGEs [11]. Note that at any given scale, 
the Ki can be rotated away through a redefinition of 
the fields L„Hi [H, H, [H]. Within the context of 
supergravity with universal SUSY breaking we find 
it natural to rotate the Ki to at Mx [33. Non- 
zero Ki are then dynamically generated through the 
RGEs, leading to naturally small neutrino masses [H, 
m. Ho, nil in agreement with observations [l^ . The 
parameters of ^e-niSUGRA are thus given by 

Mo , Ml/2 , , tan/3 , sgn(/i) , A @Mx, (8) 

where A € {Kjk, Kjk: >'ijk} : 

and we emphasise that in any given model only one 
non-zero parameter for A is chosen at Mx- In the 
no-scale mSUGRA case, although Mq = ^0 = at 
Mgut, they are non-zero at the weak scale through 
the RGE running. 

Given the additional couplings in Wy^ as well as the 
corresponding soft SUSY breaking terms, we have the 



Aijfe /^ijk ■^2jk ■^Sjk ^ijk 

0.07 0.28 0.56 0.52 pert. 



TABLE I: Weakest 2a direct bound on the five classes of 
couplings from low-energy processes such as a measure- 
ment of the Fermi constant in muon decay. Most of the 
bounds scale with (m/100) GeV, where m is the relevant 
scalar fermion mass for the process leading to the bound, 
"pert." indicates that the weakest bound only comes from 
perturbativity of the coupling Oil). Only one non-zero 
coupling at a time is considered. 



following changes in the phenomenology compared to 
Pg-mSUGRA. 

1. Through the additional couplings the RGEs 
change and thus the spectrum and couplings 
change at low-energy [H, H^] . 

2. For ^6 the LSP is unstable and the cosmological 
bound [3^ no longer applies. Thus any super- 
symmetric particle can be the LSP 

LSP e {l>i, f± ^, qL.R. X?, xf , 5} • (9) 

Within the ^g-mSUGRA model the nature of the 
LSP is determined through the initial conditions 
at Mx, Eq. and the RGEs. 

3. Through the M^^g -couplings, sparticles can be 
singly produced at colliders [48| , e.g. resonant 
slepton production at hadron colliders (49l |. 

4. Due to the possible changes in the spectrum 
as well as the additional interactions, the de- 
cay patterns of the supersymmetric particles can 
change, see for example Rcf. 50]. 

It is the purpose of this paper to investigate points 1 
and 2 in detail. Points 3 and 4 will be exemplified in 
benchmark points later, but their detailed study will 
be left for later work. First we review the experimental 
bounds on the couplings A, and consider experimen- 
tally excluded regions of parameter space. 

A. Bounds on the A Couplings 

The low-energy 2cr bounds on the tri-linear J^g- 
couplings are summarised in Refs. 0, HH, I'^ Ta- 
ble H] we show the largest allowed coupling in the five 
"classes" of couplings for scalar masses of 100 GeV. 
We do not include bounds which depend strongly on 
assumptions about quark mixing, such as those from 
RGE-induced neutrino masses, but rather those which 
come from more direct sources listed in Ref . [s^ . The 
bounds shown are to be applied at M^, and no as- 
sumption about the high energy completion of the 
model has been made. We see that the weakest A" i. 
bounds are not constraining with only perturbativ- 
ity as a limiting effect (<3.5 by some definitions). 
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d-mixing 


0.046 9.1 • 10"* 9.1 • 10"* 9.0 • 10"" 


u-mixing 


0.046 0.15" 0.15" 0.15" 



TABLE II: Weakest bound at Mx on the four classes of 
lepton number violating couplings from the RGE generated 
Ki for SPSla. The superscript x denotes cases where there 
is no neutrino mass bound on some couplings in that class 
and so we have indicated the bound coming from the lack 
of tachyons. In the case of A couplings the weakest bound 
arising from the RGEs is A231 < 0.55" from the lack of 
tachyons [53|. We have instead included the translation 
of the weak scale bound to Mx- The weak scale bound 
scales with the slepton mass (m|/100 GeV) and is stronger 
than the tachyon bounds for slepton masses below the TeV 
scale. 
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FIG. 3: Mass of ft" in no-scale mSUGRA parameter space. 
The blacked out region corresponds to points that fail neg- 
ative LEP2 Higgs constraints or that possess tachyonic 
scalars. 



It should be noted that three of the A" couphngs 
involving the first two generations have quite severe 
bounds 52]. The weakest A-^j, bounds also appear to 
aUow large ^ 0{1) couplings, but in fact many of the 
bounds are at the 10^^ level. Finally, the A^fe bounds 
are all around 5 • 10"^. 

However, within the context of a supergravity model 
formulated at the unification scale, we must take into 
account the /"e-R-GEs. Considering the effect of the 
gauge couplings on the running of an individual A 
strengthens the bounds by a factor ~ (1.5, 3.5, 4.3) at 
Mx for the couplings (A, A', A"), respectively |39l. Is^. 
The RGEs for the A's are highly coupled. In partic- 
ular, the phenomenological requirement of /i ^ to- 
gether with some A or A' ^ at My leads to non-zero 
Ki at the electroweak scale [H, HO, El] . Through mix- 
ing of the neutrinos with the neutralinos, Ki ^ im- 
plies a non-zero neutrino mass [3] . Applying the cos- 
mological bound on the neutrino mass, the resulting 
upper bounds on the Xijk{Mx), ^jki^^x) were ob- 
tained in Ref . [s^ . The corresponding weakest bounds 
at Mx on the four classes of couplings are shown in Ta- 
bleinifor the SPSla point. Here, we have updated the 
bounds using the latest SDFTSUSY version [53*]. Where 
the neutrino mass bound is not constraining, there is 
an upper bound coming from the non-generation of 
negative mass squared scalars via RGE effects. Such 
"tachyon" bounds are denoted by a superscript "x" in 
the table. Most of the bounds strongly depend on the 
origin of CKM mixing in the SM quark sector [E^, [H^l • 
We see from the table that if the CKM mixing origi- 
nates solely in the down-quark sector we obtain much 
stricter bounds than solely up-quark sector mixing. 
The nature of the bounds within the supergravity con- 
text is important for models of flavour formulated at 
the GUT scale [45)]. 



IV. BOUNDS ON fe NO-SCALE 
SUPERGRAVITY PARAMETER SPACE 

Having discussed the bounds on the A couplings we 
next consider the constraints on I^g-mSUGRA parame- 
ter space arising from precision measurements. Specif- 
ically, we consider the anomalous magnetic moment of 
the muon (17 — 2), the branching ratio Br(6 — > 57), and 
the branching ratio Br(i?s ji^ ji^). We also briefly 
discuss the impact of the LEP Higgs mass bound. 

All of the numerical calculations in the present paper 
have been performed using an unpublished Tf^ version 
of SOFTSUSY [11] ■ The predictions for the Pg conserv- 
ing Brpg(Bs fJ'~^t^~), (.9-2)^, and Brpg(6 57) are 
calculated using micrOMEGAsl . 3 . 6 [55|. Throughout 
this paper, we use the default inputs Mz = 91.1875 
GeV and nit — 172.5 GeV for the pole masses of the 
boson and top quark, respectively. The weak- 
scale gauge couplings are set to their central values 
in the WS scheme a"^(Mz) = 127.918, as{Mz) = 
0.1187. Light quark masses are also set to their 
central values in the MS scheme: mb{mb) = 4.25 
GeV, m„(2 GeV) = 0.003 GeV, md(2 GeV) = 0.00675 
GeV, ms{2 GeV) = 0.1175 GeV and mXmc) = 1.2 
GeV d^. 

The prediction for the lightest CP-even Higgs mass, 
m/i, is shown in Fig. [3] as a function of the no-scale 
mSUGRA parameter space. As discussed above, a 
lower bound of 111.4 GeV is imposed upon SDFTSUSY's 
prediction of mh. We find rah < 124 GeV for values 
of Ml/2 < 1500 GeV. 
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FIG. 4: No-scale Pe contribution to Saf^^^''^'^, i.e. Mo = 
ylo = at Mx- The blacked out region to the left and 
upper left is excluded due to LEP exclusion bounds and 
tachyons. The black contours show the region of parameter 
space that is consistent with the empirical central value of 
g — 2, plus or minus la. 

A. THE MUON (g - 2) VALUE 

The anomalous magnetic moment of positively and 
negatively charged muons is one of the very few pre- 
dictions of the SM that are not in good accordance 
with experiment. The combined experimental mea- 
surement is [13 

^oxp ^ ^ 11659208.0(6.3) x IQ-io. (10) 

The SM prediction requires the input of data on 
hadronic spectral functions. There is a well-known dis- 
crepancy between the predictions based on data from 
e+e" ^hadrons ("e"'"e~-based") and those based on 
semi-leptonic tau decays ( "tau- based" ) [58l| . Following 
the particle data group ^9|, we focus on the e+e^- 
based predictions and obtain 

Sa^ = a^^P - af^ = (22.2 ± 10.2) x lO^^", (11) 

where all theoretical and experimental errors have 
been added in quadrature. Suf^ is smaller and within 
one sigma from zero when the r-based prediction is 
used. We shall constrain the fg-MSSM such that the 
predicted 5a^^^^ = a^^P — a^^^^"*^ provides the miss- 
ing component of (g — 2)^ in Eq. (fTT|) . 

In order to compute ^a^^^^"**^, we use 
micrOMEGAsl . 3 . 6 [5a] for the one- loop super- 
symmetric Pg part, which can be enhanced by tan/?. 
The one-loop contribution is of order 

for A e {Xijk, Kjk} [s3' where m is a sparticle mass 
that appears in the loop. For m > 300 GeV and 
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FIG. 5: Pe-MSSM no-scale mSUGRA contribution to 
Br(ti S7). The blackened out region is excluded due 
to LEP measurements and the existence of tachyons. The 
contour shows the mid- value of Br(6 57). 

]A]'^ < 0.1 (as is the case for the points studied here), 

((5a^)"^6 < 0(10"^"'^) and so we neglect it, taking the 
Pg contribution only. 

Fig. (HI shows the allowed region of (5a^^^^"'^ in the 
no-scale MSSM in the Mi/2-tan/3-plane as the region 
between the contours. At higher SUSY masses (higher 
M1/2), (5a^^^'^"*'^ becomes smaller as the SUSY contri- 
butions decouple due to mass suppression of sparticle 
propagators in the one-loop diagrams. Here, the pre- 
diction of is similar to that of the SM. (5a^ increases 
for higher tan /3, which is why the contours move to the 
right with increasing tan f3. 

B. Br(6 57) 

Next we consider the measurement of Br(& — > 57). 
The current experimental result is [61] 

Br(6 ^ s7)oxp = (3.55 ± 0.26) x 10"^ (13) 

where we have added the statistical and systematic 
errors in quadrature, including the uncertainty for the 
shape function of the photon energy spectrum. If we 
now include a combined theoretical error of 0.30 x 10~^ 
[62i ] and add it in quadrature to the experimental error 
we obtain the allowed range for the central theoretical 
value at the 2(t level [1^ 

2.76 X 10"* < Br(& S7)th < 4.34 x lO""*. (14) 

The Pe-MSSM branching ratio is plotted in Fig. (O 
for the no-scale mSUGRA model. Again, A/1/2 and 
tan (3 are varied. The contour shows the mid-value of 
Br(6 -> 57), i.e. 3.55 x lO"'^. 

The Pg conserving contribution to Br(6 —>■ 57) of the 
no-scale mSUGRA model is thus in agreement with 
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the experimental bounds at 2cr in almost the entire 
plotted region. An estimate for the fg contribution 
to this process is given in Ref. [4l[. There are many 
potential contributions from products of two different 
A couplings. Although we assume only one non-zero A 
coupling at the GUT scale, several different couplings 
are generated by the RGE evolution to the weak scale 
providing a non-zero contribution. However, the 
couplings generated by the RGEs are suppressed by 
loop factors such that the I^g generated contribution 
to Br(6 — > 57) is also highly suppressed. We have 
generalised Ref. [4l| to include CKM quark mixing and 
explicitly checked that the contribution to Br (6 sj) 
is less than 10~^ in every case studied and therefore 
negligible. For example, the difference between the 
predicted values of Br(6 — > 37) in the Pg case and the 
combined Pg and J^g case for the no-scale mSUGRA 
point Ml/2 = 400 GeV, tan/3 = 20 is always smaller 
in magnitude than 10~^ for quark mixing purely in the 
down sector and for any A.^^j, set at its upper limit. A 
similar conclusion holds in the case of up-mixing. 



Br(B, 



The Tevatron experiments have recently tightened the 
upper limit on the rare decay branching ratio Br(i3s 
currently the best CDF bound is [M] 



Br, 



< 1 X 10"^ 



at the 95% C.L. The SM prediction is [H: 
Btsm{Bs /i+M") = (3-35 ± 0.32) x 10"^" , 



(15) 



(16) 



In the Pg-MSSM, the branching ratio can be signifi- 
cantly larger than in the SM since it is proportional to 
{tanPf/MX [H, [63. The tan /3 enhancement comes 



from the dependence on bottom and muon Yukawa 
couplings. The above experimental bound, Eq. (|15p . 
constrains some of the Mq 7^ mSUGRA parame- 
ter space [Hi. In Fig. ^ we show the Pg-mSUGRA 
prediction for Br(_Bs —^ fj.'^ for the specific param- 
eter range of the no-scale mSUGRA model, i.e. with 
Mq = 0, as a function of M1/2 and tan/3. We see that 
in this case, it is at least an order of magnitude below 
the experimental bound over the whole plane. 

Since the Pg-mSUGRA contribution is small, we look 
at the ^g contribution alone, i.e. without interfer- 
ence with the SM or Pg conserving contribution. The 
I*g contribution to Br{Bs — > /i+/i~) was calculated in 
Ref. [bo). Here, we generalise the calculation in order 
to include the L-R-mixings of the mediating squarks 
and also allow for non-degenerate sparticle masses. 
The computation is presented in the Appendix. We 
use this generalised result in an unpublished fg version 
of SDFTSUSY to give an estimate of the contribution to 
the decay rate. In the case of one dominant |^g oper- 
ator at the GUT scale, the contribution to Br{Bs 
/i+/i^) is always less than 10""'^^ once the A have been 
bounded by experimental constraints, irrespective of 
the nature of the CKM mixing. For our analysis, the 
fg contribution to Bt{Bs — > fj,^fj,~) is therefore ne- 
glected in the rest of the paper. However, we note that 
if we drop the requirement of only one non-zero GUT- 
scale A coupling, much larger |^g contributions are pos- 
sible. For example, choosing Xi22{Mcut) = 6.3x10"'' 
and A'j32(AfGc/T) = 7.6 x 10~^ achieves a I^g contribu- 
tion to Bi{Bs — > fj,^fi^) as high as lO"*". The corre- 
sponding product bounds for weak-scale couplings are 
given in Ref. [i^ 



V. fs-SPECTRUM 

In this section, we investigate the supersymmetric 
spectrum in the J'g-mSUGRA model, taking into ac- 
count the constraints of the previous section. We are 
particularly interested in the nature of the LSP, which 
is of fundamental importance for the resulting collider 
signatures. However, we first investigate the overall 
shift in the supersymmetric spectra due to the A cou- 
plings. Finally, we study the changes in the mass or- 
dering of sparticles. The latter determines the cascade 
decay patterns, which in turn are also essential for the 
collider signatures. 



A. Shift in Spectrum due to A-RGEs 

In order to determine the low-energy ^g supersymmet- 
ric spectrum from the parameter set at Mx, Eg. ([5]), 
we must employ the fa RGEs [H, HI, SO, lil| . They 
have been programmed at one-loop order in an unpub- 
lished version of SOFTSUSY The leading effect of 
the A couplings on the running of the supersymmetric 
masses is ~ 0[A^ /{16n^)\n{MGUT/Mz)]. Given the 
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TABLE III: Low-energy supersymmetric spectra for the 
point SPSla in the conserving case, column two, and in 
various maximum /fe violating scenarios in columns 3-5. 
In the latter other than A{Mx) all SUSY parameters are 
as in SPSla. The masses are given in GeV and we have 
high-lighted in bold face those changes compared to the 
Pe-spectrum which are 5 GeV or more. In the case of A331, 
we assume the entire CKM mixing is in the up-sector. In 
the case of A212 we have assumed the CKM mixing is in 
the down-sector. 



bounds in Tabic [TTl in most cases we only expect small 
shifts in the spectrum in the mSUGRA scenario. In 
order to illustrate the effects, we choose as a compar- 
ison the low-energy spectrum of the SPSla point [l^ 
with conserved Pg. It is listed in the second column of 
Table Hill and also shown in Fig.[T] In the first column 
we have listed the supersymmetric particles. 

In the remaining columns, we show various J'g spec- 
tra resulting from one non-zero coupling at Mx and 
with the Pg conserving parameters as for SPSla. In 
the third column of Table IIIIl we list the masses for 
the largest allowed LLE couphng: Ai23(Mx) = 0.08. 
We see that the maximal relative shift of any super- 
symmetric particle is about 2% and is in the eji mass. 
For most particles the shift is below 1%, which is in- 



distinguishable from the Pg spectrum at the level of 
significance we display. In the fourth column, we list 
the spectrum for \'^^^{Mx) = 0.122, i.e. the case of 
up- mixing [t^I- Here we obtain significant shifts in 
the masses of Vr and fi^2, which directly couple to the 
operator X^^iL^Q^Di. In particular, this leads to the 
novel possibility of a tau sneutrino LSP, which we dis- 
cuss in more detail below. There are also slight shifts 
in the di, 61, ^2 masses, with the operators again cou- 
pling directly to the dominant operator. In the 
case of down-mixing, X'^.'n is bounded to be less than 
10"'^ and the spectrum is indistinguishable from the 
Pg spectrum in the second column and we do not list 
it. 

The A" couplings are not bounded by neutrino masses 
and in some cases can still be at the perturbativity 
limit. As an example, we choose \2i2{Mx) = 0.5. We 
show the spectrum in the last column on the right. 
Here we obtain significant shifts (~ 30%) in the masses 
of Cl, di, Sl, which are all dominantly SU(2) singlets 
and thus directly couple to the leading operator. 
There are also small shifts in the corresponding dou- 
blet masses as well as the gluino mass (originating 
from changes to the 1-loop threshold contributions 
from the squarks). We expect similar effects for the 
other large A" couplings, i.e. the fields which directly 
couple should experience the largest mass shifts. 



B. Nature of the LSP 

As pointed out in Eq. in the ^g-MSSM in principle 
any supersymmetric particle can be the LSP. How- 
ever, given the restricted parameter space of Eq. ^ 
at Mjf, it is a computational question as to which 
LSPs are attained in a given model. In the general 
mSUGRA case and for vanishing A-couplings, we have 
the well-known (but usually ignored) region shown 
in the lower right of Fig. [2l where the lightest stau 
(n) is the LSP. For Aq = -100 GeV, tan/3 = 10, 
sgn(/i) = -1-1 the black line denoting the border be- 
tween the x^-LSP region, at high values of Mq and the 
fi-LSP region at low values, is approximately given by 

Ml/2 = 3.8 ■ Mo + 175 GcV . (17) 

Such a fi-LSP region always exists in mSUGRA mod- 
els, since for rising M1/2 and fixed Mq the gaugino 
masses grow more quickly than the slepton masses, 
see the approximate formulae in Ref. [l3| . However al- 
ready for Mo ~ 200 GeV, we must have M1/2 ~ 1 TeV 
for a stau LSP. We can read-off the no-scale case along 
the X-axis in Fig. [21 where Mo = 0. This always results 
in a fi-LSP. 

If we now turn on the A-couplings, this picture is in 
principle modified [i^. However as seen in the previ- 
ous section, for the maximally allowed couplings the 
shifts in the spectrum are typically modest for all but 



9 





200 






180 






160 




[GeV] 


140 




120 






100 




B 


80 






60 


staul 




40 


sneutrinoS 




20 


neutralinol 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 
coupling 
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X'33i{Mgut) ■ All quark mixing has been placed in the up 
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the A"-couplings and exceptional cases of the X'-^^^ cou- 
phng. For SPSla, the charged sleptons are the next 
lightest particles. Other than the stau, we thus might 
expect a fli LSP. However, we have checked that this 
requires Xijk- or A'-couplings far exceeding the up- 
per bounds given in Table [Til [III • (The A"-couplings 
do not significantly affect the sleptons, see Table [1111 1 
In fact, in the case of down-mixing in the weak-scale 
CKM matrix, we have found no further LSP. While 
the precise numerical bounds on the couplings are 
sfermion mass dependent and were only calculated for 
the SPSla point [ssll, the qualitative conclusion re- 
garding possible identities of the LSP holds in the re- 
gion studied in Fig. [21 

The significant shifts in the spectrum in the X212 case 
mainly affect the squarks, but they are too small to 
obtain a squark LSP. However, for up-mixing in the 
case of non-zero X'^^i{Mgut), it is possible to obtain 
a tau sneutrino LSP. As can be seen in Fig. [71 for the 
very narrow region X[}^i =^ 0.12 — 0.14, a tau sneutrino 
LSP results at SPSla. The lightest neutralino becomes 
the NNLSP for X'^^^ > 0.13. 

To summarise, we have found that there is a well- 
motivated region in the |*g no-scale supergravity model 
where the fi is the LSP. This region has hitherto been 
almost completely ignored in the literature [s^ [Tlj . 
We have also found a very special region in parameter 
space where we have a tau sneutrino LSP, i.e. we have 
one of the central results of our paper 

LSPe • (18) 

In the following, we investigate the spectrum in de- 
tail for the case of a stau or sneutrino LSP. We then 
propose several benchmark points for follow-up sys- 
tematic phenomenological analyses, including detector 
simulations. 



C. Mass Ordering in the Stau-LSP Region 

In the standard mSUGRA model, supersymmetric 
particles are produced in pairs at colliders. They then 
cascade decay to the neutralino LSP through various 
intermediate supersymmetric states. The LSP escapes 
detection, which results in the typical missing trans- 
verse energy signature for the Pg-MSSM. 

We are here interested in the mSUGRA region with 
a fi-LSP. If the A couplings are small, as indicated 
by most of the bounds, sparticles will be domina ntly 
pair-produced in the usual Pg conserving channels [T^l . 
They will then cascade decay down to the lighter su- 
persymmetric particles eventually ending at the fi- 
LSP. The nature of these cascades and thus the final 
state particles will be different from the Pg-case since 
the ordering of the spectrum can change. In particu- 
lar, the lightest neutralino will be heavier than one or 
more supersymmetric particle. Finally, the resulting 
LSPs will decay into SM particles, the flavour content 
of which depends on the dominant A coupling and 
again on the sparticle mass ordering. We thus investi- 
gate in more detail the mass ordering. 

Fig. [51 shows the f-LSP mass in the Mi/2-tan/3-plane 
for no-scale mSUGRA with sgn(/i) = +1. We shall 
mainly focus on the case of lepton-number violating 
couplings. Here in the figure, the ^g-couplings are set 
to zero, which will be a good approximation for all 
cases where they are ^0.05. The blackened out region 
to the left is excluded due to the LEP2 Higgs search 
and the absence of tachyons. In the Pg limit, LEP2 
gav e an 86 GeV lower bound on the lightest stau mass 
[56|, which is also included. 

The contours in the figure delineate regions differing 
in sparticle mass ordering. The black contour shows 
where the mass of the lightest CP-even Higgs, m/jO , is 
equal to the mass of the fi. For models where M1/2 
is smaller, i.e. to the left, rufio > mf^. The dashed 
contour shows the analogous effect for m^o = mg^ , 
i.e. to the right the neutralino is heavier than the 
lightest selectron. The neutralino is then also heavier 
than the lightest smuon, since m^^ ~ mg^ , with the 
former slightly lighter. We thus have an extensive re- 
gion, where Xi is the fourth lightest supersymmetric 
particle, the NNNLSP. This should particularly have 
an impact on slepton pair production signatures. For 
squark pair production, the heavier neutralinos and 
the two charginos which can appear in the cascade de- 
cay can in principle have significant branching ratios, 
where they completely skip the lightest neutralino di- 
rectly decaying to the light sleptons. 

It is also worth pointing out that for high M1/2 and 
high tan (3, we get the following additional changes 
compared to the SPSla spectrum. 

m^o > rrif^ (19) 

771^3 < rUf^ < ''^i>i^2 ■ (20) 
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FIG. 8: Mass of the stau-LSP in the Mi/a-tan /3-plane for 
no-scale mSUGRA. The blackened out region is excluded 
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(right). The blue contour shows where m^o 



In summary, we expect a significantly changed phe- 
nomenology in the no-scale ri-LSP parameter space, 
both compared to the Pq case but also compared to 
the hitherto extensively studied I^g case (so far usually 
called the R-parity violating case) with a neutralino 
LSP. In order to investigate this in more detail in the 
future, we propose several benchmarks which are rep- 
resentative of the possible collider signatures. 



VI. fs-BENCHMARKS 

We now turn to the description of the f^g'^SUGRA 
benchmarks. In order to fully specify the model, we 
must fix both the standard mSUGRA parameters: 
Ml/2, Mq, Aq, tan/3, sgn(/z), as well as the A cou- 
plings. We first discuss the spectra of the benchmark 
points and then go on to detail sparticle decays. 



A. Sparticle Spectra 

The first no-scale mSUGRA parameter set consists of 
a rather light sparticle spectrum: Afi/2 — 400 GeV, 
tan/? = 13 and sgn(/i) = +1. The spectrum is shown 
in Fig. [HI As can be seen by a comparison with Fig. [1] 
the spectrum is similar to SPSla. The most important 
difference is a f-LSP of mass 148.2 GeV. The selec- 
tron and smuon are the next lightest sparticles, with 
masses 160.8 and 161.3 GeV respectively. The light- 
est neutralino is the NNNLSP with a mass 161.5 GeV, 
almost degenerate with the lightest charged sleptons, 
such that the direct decays are suppressed. The spar- 
ticles are light so that there should be copious SUSY 
particle production at the LHC. We have calculated 



FIG. 9: Sparticle spectrum for no-scale mSUGRA param- 
eter set: Ml/2 = 400 GeV, tan/3 = 13, sgn(^) = -1-1 and 
A = 0. 



this spectrum in the limit of zero A and the spectrum 
in Fig. [5] will be a good approximation when all A ^ 1, 
as is the case for two of our benchmarks to be specified 
below: BCl and BC2. 
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FIG. 10: Sparticle spectrum for mSUGRA parameter set 
with a sneutrino LSP: X'33i{Mgut) = 0.122, tan/3 = 10, 
Mo = 100 GeV, Ml/2 = 250 GeV, Ao = -100 GeV, 
sgn(/i) = -1-1 and the weak-scale quark mixing is solely 
in the up sector. Numerically the spectrum is given in Ta- 
ble |TTT] We have omitted sl,cl which are almost degener- 
ate with dL,UL, respectively. We have combined ur, cr, sr 
into QR. 

For the next set of input parameters, we pick a point 
from Fig. [7| with an SPSla-like spectrum but with a 
tau sneutrino LSP, in order to obtain different sig- 
natures, i.e. here we must choose a non-zero A: 
A331 (Mgc/t) = 0.122, tan/3 = 10, Mq = 100 GeV, 
Ml/2 = 250 GeV, Aq = -100 GeV, sgn(/i) = +1 and 
the weak-scale quark mixing is solely in the up sector. 
The spectrum is displayed in Fig. 1101 Numerically it 
is given in Table IIIIl The mass ordering of the light- 
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est sparticles is nif-^ > m^o > mp^ . The spectrum is 
light and would lead to copious SUSY production at 
the LHC. 

The final no-scale mSUGRA parameter set consists 
of a somewhat heavier spectrum: Afi/2 = 600 GeV, 
tan/3 = 30, and sgn(/Lt) = +1, X212{Mgut) = 0.5. 
The spectrum is displayed in Fig. [11] note the scale 
difference to Fig. M "SUSY" detection and measure- 
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FIG. 11: Sparticle spectrum for no-scale mSUGRA pa- 
rameter set: Ml/2 = 600 GeV, tan/? — 30, sgn(/i) = +1, 
Ai.'i2(MGC/T) = 0.5. 

ment at this point will be more difficult than for the 
no-scale set I, but still easily possible Q. We have 
chosen this point to represent the novel mass ordering 
where the neutralino is the NNNLSP, as discussed in 
the previous section: m-^o > mg^.^^ > rrif^. 



B. Sparticle Decays 

As pointed out in Ref . [s^ , depending upon the flavour 
structure of the A couplings, the non-neutralino LSP 
may dominantly undergo 2 or 4-body decays. Since 
we expect two LSPs per SUSY event, the 4-body de- 
cays will significantly change the number of particles 
in the final state and therefore alter the potential sig- 
natures. A benchmark where 2-body LSP decays dom- 
inate and another where 4-body LSP decays dominate 
then seems expedient. We choose a stau LSP for the 
4-body decay, since the calculations for the partial 
widths have been completed W\ and the sneutrino 
has no 4-body decays. We choose a point with a tau 
sneutrino LSP, which necessarily has 2-body decays 
because of the dominant A33]^ coupling, which it di- 
rectly couples to. 

For small A, the LSP may be long-lived, leading to 
sparticles with a measurable detached decay vertex in 
the detector. If the |g couplings are even smaller, such 
decays are longer lived than ~ 10~^ seconds, the staus 
will decay outside of the detector and appear as heav- 



ily ionising tracks. This possibility has been examined 
recently in the literature [ssj , albeit in a different con- 
text. In an intermediate coupling regime, signifi- 
cant numbers of the staus will decay in the detector 
providing heavily ionising tracks with detached ver- 
tices. This possibility is new and interesting and so 
we shall include a benchmark that covers it. 

We note in passing that the case of a long-lived LSP 
in principle also include timescales which are relevant 
for cosmology, i.e. t(LSP) > O(lsec) and the result- 
ing bounds must be taken into account [755. As we 
saw, in the case of R-parity violation, in principle any 
supersymmetric particle can be the LSP. We found in 
our restricted mSUGRA scenarios either a neutralino, 
a sneutrino or a stau LSP. As has been recently in- 
vestigated in some detail the cosmological bounds on 
a long-lived charged particle for example from nucle- 
osynthesis are very different for charged 7^ 82 ^ sl, ^ 
and neutral particles [zl, [li H H IZl, M l8l|, |8l|. 
For example the long-lived charged particles can form 
bound states with the free nuclei thus affecting their 
nuclear reactions, which in turn directly affects the 
abundance of the light elements 83] . A more detailed 
analysis is unfortunately beyond the scope of this pa- 
per. 

^ijk couplings necessarily lead to 4-body stau decays. 
These involve many jets in the final state as well as 
taus from the f decays. Thus the fourth benchmark 
shall have a non-zero A^^^, at Mqut- More precisely, 
the four benchmarks chosen are 



• BCl: no-scale mSUGRA with \i2i{Mgut) = 
0.032, tan/3 = 13, M1/2 = 400 GeV, sgn(^) = 
-|-1. The spectrum is shown in Fig. [51 

• BC2: no-scale mSUGRA with Wr^-^-^^Maur) = 
3.5 X lO-'^, tan/3 = 13, M1/2 = 400 GeV, 
sgn(/i) = -|-1. The spectrum is shown in Fig. [9l 
The LSP is long-lived. 

• BC3: mSUGRA with \'^^^{Mgut) = 0.122, 
tan/3 = 10, Afo = 100 GeV, M1/2 = 250 GeV, 

= —100 GeV, sgn(^) = -1-1 and weak-scale 
quark mixing solely in the up sector. The spec- 
trum is shown in Fig. [TUj and given in Table IIIII 

• BC4: no-scale mSUGRA with A212 (^GC/t) = 
0.5, tan/3 = 30, M1/2 = 600 GeV, sgn(/i) = +1. 
The spectrum is shown in Fig. 1111 

In order to calculate the decay rates of sparticles, 
we pipe the output of the SDFTSUSY spectrum cal- 
culation through ISAWIG1.200. This is linked to 
ISAJET7.64 [8y| in order to calculate the 2-body par- 
tial widths of the sparticles and Higgs'. The output 
from this procedure is fed into a specially modified 
version of the HERWIG program [s^l that is able to cal- 
culate partial widths for the 2- and 4-body stau decays 
[111. We shall now display the decay branching ratios 
for sparticles and Higgs'. We omit decays since 
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these follow those of the Standard Model limit very 
closely. 

The decays for benchmark point BCl can be seen in 
Table IIVI Here and in the following tables of bench- 
mark points, for decays into the light quark flavours 
the results are summed over the flavours q = u,s, d, c. 
We see that, as designed, the fi-LSP is completely 
dominated by 4-body decays. Each 4-body decay is 
accompanied by a tau. Thus good tau identification 
would help to identify this scenario. Furthermore, 
each 4-body decay includes a final-state neutrino, re- 
sulting in missing transverse momentum, ^t- This 
should however be reduced compared to a Pg model, 
since the stau momentum is diluted in the 4-body de- 
cay. Also, the left-handed selectron and sneutrinos 
undergo significant decays. The heavier neutrali- 
nos have significant branching ratios into charginos as 
well as Z'^-bosons/Higgs and lighter neutralinos. The 
number of expected taus from a given SUSY pair pro- 
duction event is often 4 or more. We also see a differ- 
ence in the decays of and [Ir due to the presence of 
the A121 coupling. Furthermore, we see from the table 
that other J*g couplings that are induced in the RGE 
running are not large enough to induce branching ra- 
tios greater than 0.005. 

The decays for benchmark point BC2 can be seen in 
Table IVl We see that, as designed, the fi-LSP is com- 
pletely dominated by 2-body decays into non-& jets. 
Unlike BCl, the LSP decay comes without a neu- 
trino, meaning that SUSY events do not necessarily 
possess the classic missing transverse momentum 
signature. A comparison between Tables |Vl and IIVI 
shows that many of the decays are identical as a result 
of the spectrum being approximately identical. The 
main differences are in the light sparticle decays: fi, 
^L,R, P'L,R and Ve^fj,. The X'^n coupling is too small to 
have a significant affect upon squark decays. However, 
in a collider the main difference will be the existence 
of detached vertices from the relatively long-lived fi, 
which has a lifetime of Tf ^ = 1.1 x 10""'^^ sees in its rest- 
frame. This corresponds to a rest-frame decay length 
of CTfj w 0.3 mm. However, when this is boosted into 
the lab frame (for example by a factor 7 — 30), then 
the decay-length becomes about 1 cm. If these f-LSPs 
come from other sparticle decays, once can expect a 
coincidence of detached vertices with other particles 
originating from the interaction point in the detector. 
Of course the number of decays will fall approximately 
exponentially with radial distance, and so one may ob- 
tain decays at radial distances of /xm-mm, possibly in- 
terfering with 6-tagging. We leave any in-depth study 
of such effects to future work using the benchmarks. 

The decays for benchmark point BC3 are shown in Ta- 
ble lVIl Here the tau sncutrino is the LSP and it neces- 
sarily couples to the dominant operator. Therefore 
it will always decay via the 2-body mode to the purely 
hadronic final state bd (bd). One of the daughter jets 
involves a bottom quark, providing for the possibil- 



ity of using 6— tagging in order to help identify SUSY 
events. The neutralino is the NLSP and the only de- 
cay mode open is the 2 body decay to Vt-Dt {i>*Vt), 
which results in missing in the final state. Both 
xf and X2 have significant branching ratios to tau lep- 
tons/stau sleptons. Thus the cascade decays of left- 
handed squarks and indirectly thus also the gluinos 
should still lead to significant number of tau's in the 
final state. We also see from the table that other Pg 
couplings that are induced in the RGE running are 
not sufficiently large to induce branching ratios greater 
than 0.005. 

We display the decays for benchmark point BC4 in 
Table IVIII Here, the fi-LSP decays exclusively into 
a r and 3 non-6 jets, since the decays via a virtual 
chargino are strongly suppressed, since the tau and 
the daughter quarks are dominantly SU{2) singlets. 
Thus, this point will provide very little signal jJt- The 
first two generations of right-handed squarks undergo 
dominant J*g decays into two jets, altering many LHC 
signatures. Aside from these, the decays are rather 
similar to those of BC3, as a comparison with Ta- 
ble I VII shows, although some quantitative differences 
are present in every channel. The large number of jets 
in final states from this model could make analysis of 
SUSY events difficult: combinatoric backgrounds are 
likely in event reconstruction. Examining specific de- 
cay chains involving Icptons can help reduce these (soj . 
The combinatoric backgrounds problem will become 
worse for high luminosity running, where pile-up will 
increase the number of jets in each event. However, 
resonant squark production ought to be possible, pro- 
viding an additional empirical handle on the model. 



VII. CONCLUSIONS 

We have investigated the spectrum of the I*g-MSSM 
embedded in supergravity, including indirect con- 
straints on parameter space. We have found different 
regions of parameter space with a neutralino, a stau 
and also a tau sncutrino LSP. Taking these into ac- 
count, we have presented the first set of benchmarks 
in I*g mSUGRA. All of the points are within current 
search limits and are consistent with measurements 
of precision observables. We have picked different 
SUSY breaking scenarios: light and heavy. The heav- 
ier benchmarks are more difficult to detect at the LHC 
and so should be used to see how much is possible 
to achieve through data analyses in this difficult sce- 
nario. The light benchmarks are designed to enable 
high statistics analyses in order to determine what is 
possible to divine from LHC data in the optimistic 
case. 

We have picked the fiavour of the J*g couplings in order 
to show-case various features relevant for the experi- 
ments. In the first benchmark, four-body decays of the 
LSP are expected including neutrinos, leading to com- 
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TABLE IV: Sparticle and Higgs decays for BCl: A121 (Mgu 
sgn(/i) = +1. |*g decays are shown in bold font. Only decays 
decays are prompt. The spectrum is displayed in Fig. (9) 
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) = 0.032, tan/3 = 13, M1/2 = 400 GeV, Mo = Ao = 0, 
nth branching ratios greater or equal 0.5% are shown. All 



plicated final states and ^t- The second benchmark 
has vertices detached from the interaction point, where 
LSPs decay into non-bottom jets. Coincidence with 
SM particles from higher up the cascades are expected 
from the interaction region. In the third, 2-body de- 
cays of a tau sneutrino LSP into jets (one of them a 
bottom jet) are expected. The fom'th benchmark has 
LSPs promptly decaying into 3 jets and a tau. It also 



has a large coupling capable of producing significant 
single-squarks at hadron colliders such as the LHC. 
In order to enable analysis of the proposed bench- 
mark points, we have provided HERWIG files on the web 
at URL address http://hepforge.cedar.ac.Uk/~ 
allcLnach/rpv/, where they are available for public 
download. 
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TABLE V: Sparticle and Higgs decays for BC2: \3ii{Mgut) = 3.5 x 10"^ tan/3 = 13, Mj/j = 400 GeV, Mo = Ao = 0. 

decays are shown in bold font. Only decays with branching ratios greater than 0.005 are shown. In the tenth row 
qR = d-R, sr, ur, cr. And the q in the decays is correspondingly q — d,s, u, c. All decays are prompt except for fi, which 
has a hfetime of 1.1 xlO"^^ sec. The spectrum is shown in Fig. (5] 
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TABLE VI: Sparticle and Higgs decays for BC3: A331 (Mge/t) = 0.122, tan/3 = 10, M1/2 = 250 GeV, Mo = 100 GeV, 
^0 = —100 GeV. 1*5 decays are shown in bold font. Only decays with branching ratios greater than 0.005 are shown. All 
decays are prompt. 
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TABLE VII: Sparticle and Higgs decays for BC4: A212 (Mgc/t) = 0.5, tan/3 = 30, M1/2 = 600 GeV, Mo = ylo = 0. f, 
decays are shown in bold font. Only decays with branching ratios greater than 0.005 are shown. All decays are prompt. 
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APPENDIX A: CONTRIBUTION TO 



The final expression for the partial width of the decay 



B. 



e+Cj reads: 



Here, we generalise the calculation of Ref. [69| to 
include non-degenerate sparticle masses as well as 
squark mixing. The diagrams of interest consist of 
two s-channel diagrams mediated by sneutrinos, and 
a t-channel diagram mediated by left handed up-type 
squarks. Using the convention of [53i] . we write the 
squark mixing as: 



ULj = CjUlj — s 



J 



(Al) 
1, 2, 3 is a 



where the Cj = cosOJ, sj = sin0" 
generation index and the states on the right hand side 
are mass eigenstates. For convenience, we define the 
following quantities: 



Im 



3 \ * \/ 

ml 

3 \ * \ / 

ml 

1 ^ 2 

2 ^li]^mi3[^2 



(A2) 

(A3) 

(A4) 
(A5) 



■^Bq- ^Bq. 

x{[(|ALP + |5LP)(i--?--L 

2Re(A^„,C^,*)Ac„(l + «2_^^) + 



(A6) 



where . , Mb^ . are the decay constant and mass of 
the Bq. meson respectively. 
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